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ABSTRACT 



Context. Following the recent insight in the material structure of comets, protoplanetesimals are assumed to have low densities and 
to be highly porous agglomerates. It is still unclear if planetesimals can be formed from these objects by coUisional growth. 
Aims. Therefore, it is important to study numerically the collisional outcome from low velocity impacts of equal sized porous ag- 
glomerates which are too large to be examined in a laboratory experiment. 

Mefhods. We use the Lagrangian particle method Smooth Particle Hydrodynamics to solve the equations that describe the dynamics 
of elastic and plastic bodies. Additionally, to account for the influence of porosity, we follow a previous developed equation of state 
and certain relations between the material strength and the relative density. 

Results. Collisional growth seems possible for rather low collision velocities and particular material strengths. The remnants of colli- 
sions with impact parameters that are larger than 50 % of the radius of the colliding objects tend to rotate. For small impact parameters, 
the colliding objects are effectively slowed down without a prominent compaction of the porous structure, which probably increases 
the possibility for growth. The protoplanetesimals, however, do not stick together for the most part of the employed material strengths. 
Conclusions. An important issue in subsequent studies has to be the influence of rotation to collisional growth. Moreover, for realistic 
simulations of protoplanetesimals it is crucial to know the correct material parameters in more detail. 

Key words, planetary systems; formation - planetary systems: protoplanetary discs 



1. Introduction 

The initial growth of particles in a protoplanetary disc is accom- 
plished by sticking collisions. All solid objects in a planetary 
system are believed to have been developed from collisions be- 
tween small dust particles with initial sizes of about 0.1 micron 
(iGreenberdI 19821) . These collisions result from Brownian mo- 
tion, gasjurbu^^ 

fT977t IVoelk eTaniT980l: IWeidenschilling & Cuzziiri993l) . It has 



■ been shown that even protoplanetesimals of equal mass might 
have rather high relative velociti es due to different shapes and 
hence different gas resistances (Benz 2000), depending on the 
orbital distance to the protostar. Experimental data for colli- 

, sional growth of macroscopic agglomerates with sizes extend- 
ing the dm-regime are rather scarce. The formation of such 
macrosco pic bodies from s maller dust grains was studied in 
detail by Ossenkopf" (1993) a nd the frac tal growth w as mod - 
elled by Dominik & TielensI (Il997h and iKempf et all (Il999h . 
Generally, the dust particles stick together due to van-der-Waals 
forces. This growth process leads to (highly) porous agglomer- 
ates^ fact that was experimentally confirmed by Wurm & Blui3 
(119981) . These macroscopic agglomerates do not stick together 
as easily as tiny dust grains which only interact by van-der- 
Waals forces. Instead, aggregate compaction, fragmentation and 
disruption become important above a specific kinetic energy 
of the collision. The agglomerates can break and can be dis- 
persed which eventually prevents a fast growth process or 
even further growth at all. It is now generally assumed that 
decimetre sized porous agglomerates can form rather quickly 
in th e presolar nebula and the protoplanetary accretion disc 
JPopp e et al. 2000; Wurm & Bl um. 1998; Blum & Wurm .2000; 
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IWeidenschiiiingl2000h . The critical bottleneck in the formation 
process of planetesimals lies however in the size range from 
about several decimetre to objects whose interaction is mainly 
dominated by gravitation, that is 1 to 10 km planetesimals. First 
studies to the collisional growth of brittle planetesimals with ma- 
terial stren gth we re performed by Benz & Asphaug ( 1994, 19951 
19991) and iBenzl (|2P00). Their results indicate that the weakest 
bodies during collisions are in the size range from metres to 
around 100 m. Caused by the stronger interaction with the gas, 
metre sized bodies experience higher collision velocities and are 
also rather fragile objects and easily disrupted. 

A better understood problem is the formation of planets from 
planetesimals. The formation of planetary embryos or co res has 
been i nvestigated in a series of publicati ons (Stewart & Wetherilll 
119881; IWetherill & Stewart! 1 19891 Il993h with the application of 
a statistical approach to describe the evolution of a swarm 
of planetesimals around the star Later, by the use of mod- 
ern supercomputers, direct n-body simulations of a large num- 
ber o f planetesimals became fea sible and yielded similar re- 
sults dKokubo & Idal l2000l l2002h . Once planetary cores have 
formed, the terrestrial planets grow through collisio ns be- 
tween them (.C hambers & Wetherill 1998; Agnoret al. 1 19991; 
iKominami & Idall2002l) . In this size regime, the collisional out- 
come is by far domina ted by gravitational interactions (e.g., 
lAgnor & Asphaugll2004l) . 

One alternative to avoid the bottleneck at the metre length 
scale of the formation process of planetesimals is gravitational 
instability of the dust layer in the midplane of the protoplane- 
tary accretion disc. This idea was suggested independently by 
Safronov (1969) and Gold reich & Ward (1973). For the solar 
nebula, iOoIdreich & Ward ( il973i) found gravitational clumping 
of the dust which forms clumps with radii of 10 km at 1 AU 
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within 1000 years after the onset of the vertical settling of dust. 
However, it is still unclear and strongly discussed in the commu- 
nity whether the dust in the midplane can become gravitational 
unstable, because Kelvin-Helmholtz instability and turbulent 



motion of dust particles in the disc lead to significant velocit 
dispersions which may prevent the collapse (fcuzzi et al. 19931; 
Weidenschilling||200a lYoudin & Shul l200l lYoudin & Chiand 



I 



20041) . 



Therefore, collisional growth of metre sized bodies has to be 
studied in greater detail in order to investigate the mechanisms 
of planetesimal formation. Labo ratory experiments o f collisions 
between water ice objects (e.g.. [Bridges et aLlll996) have been 
restricted to smaller sizes (up to several cm) and low impact 
velocities (up to several cm/s). Thus, experiments have con- 
centrated on the dependence of sticking forces from the frost- 
coating of the surfaces (Supulveret al. 1997) or on measuring 
restitution coefficients dSupulver et al.lll995b . 

Since real collisions between objects of metre size are not 
feasible yet in laboratory experiments, one has to resort to nu- 
merical simulations. Here, we employ the numerical Lagrangian 
particle method Smoot h Part i cle Hy drod ynamics (SPH) which 
was fir st introduced by iLucvl (Il977h and iGingold & MonaghanI 
(1 19771) for modelling compressible flows in astrophysical prob- 
lems. The method was extended to solid state mechanics in 
the beginning of the nineties by iLibersky & Petschek (1990) 
and irnproved extensively in the following years (Libersky et al.. 
1991 iRandles & Liberskvl \l99^ iLiberskv etal l Il997h . 



After 



Bend (12000') has applied the method successfully to simulations 



of low velocity collisions between brittle objects, Sirono (200|) 
used experimentally measured data to parametrise the compres- 
sive strength curve of icy grain agglomerates and derived a mod- 
ified Murnaghan equation of state (see sect.|2|i to simulate porous 
objects. 

Rece nt astronom i cal inv estigations such as the Deep Impact 
mission ( Lisse et aTl l2006l) and observations of comets (e.g., 
iLamy et al.ll2006l) yield strong indications for rather low densi- 
ties of comets (the values range from 0.1 g/cm^ to several tenths 
of g/cm^ and still have some larger errors). This matches fairly 
well with the idea of collisional growth of fractal porous dust 
agglomerates to planetesimals. 

In order to extend the investigations bvl Sironol (|2004|) . we ap- 
ply his model to simulate the collisions between equal sized (ice) 
agglomerates with different impact parameters and also explore 
the influence of material properties to the simulation outcome. 
The outline of the paper is as follows: In the next section, we 
will present the physical model, that is the basic equations and 
material properties we assume. In sect. [3] we will discuss some 
numerical issues and in sect.|4]the setup of our simulations. We 
present our results in sect.|5]and finally draw some conclusions. 



2. Physical Model 

2.1. Basic Equations 

The system of partial differential equations that describe the dy- 
namics of an elastic solid body is given by three equations. The 
first is the continuity equation 



dp dv" 
df ^dx" 



(1) 



where g denotes the density and v the velocity and where the 
Einstein summation rule is applied. Greek indices denote the 



spatial coordinates and run from 1 to 3. The second equation 
in the system accounts for the conservation of momentum 



dv" 

IT 



1 do-"!^ 
g dxP 



(2) 



The stress tensor cr is given by the pressure p and the deviatoric 
stress tensor S according to 



(3) 



In contrast to fluid dynamics, this set of partial differential 
conservation equations is not sufficient to describe the elastic 
body, since the time evolution of the deviatoric stress tensor is 
not yet specified. The missing relations are called the constitu- 
tive equations which describe the rheology of the body and relate 
the kinematic states of the body to the dynamic states. The elas- 
tic behaviour of a solid body can be described by Hooke's law, 
which reads in three dimensions 



5"^^ ~2yLiJe"^- ^^^^e^'^j. 



(4) 



where p is the shear modulus of the material, and e"^ denotes 
the strain tensor which is given by 



1 Idx" 



dx'^ 



2 \ dxP dx" 



(5) 



Here, the primed coordinates denote the coordinates of the de- 
formed body. 

The stress rate has to be defined in a way that obeys the 
principle of frame invariance. There are various possibilities to 
achieve this . We follow the usua l approach which is used for 
SPH codes d Benz & Asphaug|ll994i) and adopt the Jaumann rate 
form, where the time evolution of the deviatoric stress tensor can 
be expressed as 



1 



dS"/^ 
df 

where R"^ denotes the rotation rate tensor 



R 



2 \ dxP dx" 



and the strain rate tensor 



dv" dv^ \ 
2 ^ 'd^]' 



(6) 



(7) 



(8) 



The closure of this set of equations is provided by the equa- 
tion of state that relates the pressure p to the density of the 
agglomerate. Here, we focus on porous o bjects and follow the 
semi-empirical approach by ISironol (12004 ). We use an exten- 
sion of the Murnaghan equation of state which accounts for the 
change of porosity and reads 



, , f K(e'Q)(g/g'Q -1), g< g^. 



(9) 



The quantity is the reference density of an agglomerate at no 
external stress while go is the initial density of the porous ag- 
glomerate. Note that g'^ is in general different to the initial den- 
sity go of the porous agglomerate after the start of the simulation. 
The critical density is a function of the reference density g'^ 
and determines the state where the pressure reaches the compres- 
sive strength limit l.{g). As soon as the pressure decreases, the 
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behaviour of the agglomerate is again elastic but with a differ- 
ent bulk modulus Kig'^^) since the reference density has changed. 
The slope of Kig'^^) has to be known either by experimental data 
or theoretical considerations. The tensile strength of the mate- 
rial is determined accordingly and limits the tension for negative 
pressure. 

The sound speed of the agglomerate is given by the bulk 
modulus of the agglomerate and the reference density accord- 
ing to the relation 



c(q'o) = ^K(g'^)/g'Q. 



(10) 



For plastic states when the pressure exceeds the compressive 
strength of the material the sound speed is calculated according 
to 



c(g) 



dX(e) 
V dg ■ 



(11) 



ISironol (|2004|) additionally uses an empirical damage model 
for porous agglomerates which is based on the damage mo del for 
brittle materials developed for SPH by Benz & Asphaug ( 11994 
[19951). However, this model is only apphcable for the simula- 
tion of brittle fracture in rocks dGradv & KippI [198(1) . and our 
test simulations show that the application of the damage model 
to porous agglomerates does not yield reliable results since the 
model also includes compressive damage effects. Because of 
these considerations, we do not include any damage model in our 
calculations and only consider the fragmentation due to plastic 
flow. 



2.2. Material properties 

The shape of the compressive and tensile strength curves have 
to be known either by experimental data or theoretical c onsid- 
eration s. Since we want to compare our calculations to ISironol 
(l2004l) . we primarily use the values published in his paper 
which are derived from e xperimental data (see iKendall et alj 
119871: IValverde et al.|[T998[) : The compressive strength is given 
by l.{(f) - Soi/j", with tp - g/go, and the tensile strength ac- 
cordingly by T((p) = Totfi/^. The dependence of the bulk modulus 
of the aggregate on the porosity reads K{ip) = /Tqi^''. The shear 
modulus is assumed to be yu - K/2, and the shear strength is 
defined by Y(ip) - ^j2T{tp)'Yj{ip)l3. The parameter ip is related to 
the porosity = gjg^ in the following way 



Q_ 

go 



00 ' 



(12) 



where g^ denotes the density of the solid material in the porous 
body and po is the porosity of the body at the initial simulation 
time. 

The transition from an elastic state to a yield state of a mate- 
rial can be characterised by a complex stress state. The stress at 
which the transition starts is called the yield stress and the con- 
dition is called the yield criterion or the plas ticity criterion. In 
this pa per, we follow the implementation by i Benz & Asphaud 
( Il995h and assume isotropic material and apply the von Mises 
yield criterion: We calculate the second irreducible invariant of 
the deviatoric stress tensor J2 - SapS"^ and reduce the devia- 
toric stress when necessary according to 



S-l^^fS-l^, /^min [7(^)2/372,1]. 



(13) 



3. Numerical Issues 

Smooth Particle Hydrodynamics (SPH) is perfectly suitable for 
the simulation of brittle and plastic materials. The continuum 
of the solid body is discretised into mass packages which are 
called particles. The particles move like point masses according 
to the Lagrangian form of the equation of motion. They carry 
all physical properties like mass, momentum and energy of the 
part of the solid body they represent. The particles interact by 
kernel interpolation during the simulation to exchange momen- 
tum and energy. For a complete description of the method and 
its featu res an d qualitie s, we refer to t wo comprehensive review 
articles dBen^ 1990; M onaghanll992h . 

In standard SPH, the velocity derivatives in eqs. ^ and (O 
for the determination of the rotation rate and the strain rate tensor 
for particle / are usually calculated according to 



5/ ^ ej ' 



5/' 



(14) 



where the sum runs over all interaction partners j of particle i, 
and W'^ denotes the kernel for the particular interaction. 

This approach, however, leads to erroneous results and does 
not conserve angular momentum due to the discretisation error 
by particle disorder in simulations of solid bodies. This error can 
be avoided by constructing and applying a correction tensor C^^ 
according to 



dx] 



where the correction tensor C^^ is the inverse of 



that is 



(15) 



(16) 



(17) 



By applying this correction tensor first order consistency can be 
constructed where the errors due to particle disorder cancel out 
and the conservation of a ngular momentum is ensured. This is 
similar to an approach that lBonet & Lokl d 19991) proposed for the 
conservation of angular momentum, where all spatial derivatives 
are corrected according to eq. ( fTSl l. We have found however that 
it is sufficient to correct only the rotation rate and the strain rate 
tensor 

We additionally use the standard SPH artificial viscosity 
dMonaghan & Gingold|[T983l) with a - \. However, in order to 
reduce the viscous energy dissipation, we follow Sirono's ap- 
proach and apply the artificial viscosity only when the approach- 
ing relative velocity of two particles is larger than the sound 
speed in contrast to the standard case where the artificial vis- 
cosity acts always on approaching particles^ 

Our SPH code parasph (Hipp & Ros enstielll2004l) has been 
successfully tested with some numeri cal standard tests such as 
the co llision of perfectly elastic rings dSwegldl 19921: iMonaghanI 
I2OOO ) and the simulations of ductile and britt le rods under ten- 
sion (iBenz & Asphaug|1994l:lGrav et al.l200Th . Additionally, the 
Siron o model in the extended code has been checked carefully 
(Schafer 2005) and we have s uccessfully rep roduced the impact 
result for the standard setup in lSironol d2004t) . 
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4. Setup 

All collisions have been simulated in three dimensions. The 
setup of our calculations includes the following parameters: 
Two spherical ice agglomerates with the initial density go - 
0.1 g/cnr' and identical radius r = 1 m collide at either 20 m/s 
or 10 m/s. The initial density implies an initial porosity of 90 %. 
The density dependenci es of compre ssive and tensile strengths 
are chosen according to lSironol (|2004) to l,(g) - 'Lo(g/go)^ and 
T(g) - T()(g/go)^, and the bulk modulus is given by K{g) - 
K(){g/go)'*. We always assume the shear modulus is given hy fi - 
K/2. For the first simulations with impact velocity 20 m/s we 
use the basic setup values from Sirono and set Kq - 6 x 10^ Pa, 
So - 600 Pa, and To = 6x10^ Pa. These are the values for which 
Sirono (and we) find perfect sticking behaviour for a smaller 
sphere impacting into a larger object (with a radii ratio of 10:3) 
at speed 10 m/s. The temperature implicitely assumed in this set 
of parameters coiTesponds to very cold ice, which is consistent 
with the simulation results where large deformations but hardly 
any compression is found. Since the probability of a head-on 
colUsion is low, we have additionally varied the impact parame- 
ter of the collision, b = 0.2, 0.4, 0.8, 1.2 m. The results of these 
simulations are shown in sect. 15.11 

To study the importance of realistic material properties, we 
have then varied the tensile and compressive strength parameters 
To and Eq in some of the 10 m/s head-on collisions and present 
these calculations in sect. 15.21 Reducing the impact velocity also 
allows to compare the influence of the latter. 



5. Results 

5.1. Varying impact parameter 

As noted above, in this section w e focus on the material pa- 
rameters for which ISironol (12004 ) found perfect sticking for a 
collision with a projectile 3:10 smaller in radius than the target. 
Figs-HElpresent the simulation outcome for five different impact 
parameters. The plots show the SPH particles with a colour-scale 
cod43 for the density of each particle. As we are mainly inter- 
ested in the final configuration of the impact after compression 
has ceased, the density is acting as a measure for any permanent 
compaction that may have taken place during the compression 
phase. 

Each sphere consists of approximately 1 1 500 SPH particles 
which are initially placed on a tetrahedral grid to maximise the 
number of adjacent and interacting particles and thus to increase 
the resolution. The figures show the surface of the resulting ob- 
jects after impact, as the particles are plotted opaque. 

The first simulation (see fig. [1]) demonstrates the behaviour 
during a head-on collision with impact par ameter b = . In con- 
trast to the sticking mechanism found by Sironol (12004') with the 
same material properties but different sizes, the two spheres do 
not penetrate. In fact, they undergo a large change in shape and 
the two spheres end up as flattened discs. Almost the entire ki- 
netic energy of the impact is dissipated by plastic deformation, 
and ultimately the two discs are at rest in the barycentric sys- 
tem. Only some particles which were tossed outwards during 
the compression phase fly away at constant speed. 

The picture slightly changes for the simulation with b - 
0.2 m (see fig. |2|. As in the head-on collision, the two spheres 



' As in the printed version of this paper the figures have to be dis- 
played in grey-scale, we refer the reader to the online- version for colour 
figures. 




* 
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Fig. 1. Head-on collision with impact velocity 20 m/s. The ma- 
terial strength coefficients are chosen to So = 600 Pa and 
Tq - 6000 Pa. The plot shows colour-scaled the density for each 
particle at simulation time f = 0.662 s after contact. The colour- 
bar is in units of g/cm^ and is valid for all colour-scale plots in 
this paper. Also note that the line of approach is tilted by ~ 30° 
from the vertical in all plots. 




Fig. 2. Collision with impact parameter b - 0.2 m. The plot 
shows colour-scaled the density for each particle at simulation 
time t - 0.662 s after contact. 



transform into flattened discs with some separate particles float- 
ing away. In contrast to the head-on collision, the remnants of 
the spheres rotate slowly. 

First major difference arises at the step from b - 0.2 m to 
b = 0.4 m (see fig. |3]l, when the spheres break up. Again, the 
spheres undergo a large plastic deformation during the compres- 
sion phase and become strongly prolate. A small piece is quar- 
ried out of each object after the compression due to the strong 
tensile forces caused by rotation. 

Even more smaller fragments emerge after the collision with 
impact parameter b - 0.8 m (see fig. |4]l. However, now the 
largest remnants of the collision consist of two rotating half- 
spheres. The density of the half-spheres is exactly the density 
of the initial spheres. Essentially they have not been compacted. 

The situation hardly changes for b = 1.2 m (see fig. |5]l- 
Again, smaller fragments are tossed out from the spheres af- 
ter the compaction phase. The two formed larger remnants have 
again the shape of half-spheres, slightly larger than in the simu- 
lation with b - 0.8 m. The rotation rate and the velocity of the 
half-spheres is larger than for smaller impact parameters. 
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Fig. 3. Collision with impact parameter b - 0.4 m. The plot 
shows colour-scaled the density for each particle at simulation 
time t = 0.662 s after contact. 




Fig. 4. Collision with impact parameter b - 0.8 m. The plot 
shows colour-scaled the density for each particle at simulation 
time t = 0.662 s after contact. 




Fig. 5. Collision with impact parameter b - 1.2 m. The plot 
shows colour-scaled the density for each particle at simulation 
time t = 0.662 s after contact. 
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Fig. 6. Ratio Q between kinetic energy at end of the simulation 
and initial energy for five different impact parameters. 




Fig. 7. Head-on collision with impact velocity 10 m/s. The ma- 
terial strength coefficients are So = 600 Pa and Tq = 6 000 Pa. 
The plot shows colour-scaled the density for each particle at sim- 
ulation time t = 0.5 s after contact. 



The remnants (except for a few particles) of the collision are at 
rest. Although they do not stick together, their relative velocity is 
zero. Even for impact parameters b - 0.2 and 0.4 m, the kinetic 
energy at the end is less than 10 % of the initial kinetic energy. 
The spheres are effectively slowed down during the compression 
phase. 

The situation changes for higher impact parameters. For b = 
0.8 m less than 80 % and for b = 1.2 m even less than 50 % of 
the initial kinetic energy is lost due to the compression. 

Interestingly, the density of the remnants after the simula- 
tion is basically the initial density for all studied impact parame- 
ters. Only some small areas in the contact region have a slightly 
higher density. The density in the interiour of the objects does 
not differ distinctly from the density at the surfaces. Therefore 
only the latter is shown in the figures. 



The strong influence of the impact parameter b on the change 
in kinetic energy with the respect to the centre of mass is shown 
in fig. m where the ratio Q of the kinetic energy at the end of 
the simulation to the initial kinetic energy is plotted. In the case 
of the head-on collision, more than 99.7 % of the initial kinetic 
energy has been dissipated during the plastic deformation phase. 



5.2. Varying material strengtli 

Although lSironol (l 2004l) has investigated the influence of the ma- 
terial strength on the collisional outcome, we have additionally 
studied the behaviour for our equal-sized agglomerates. 

All simulations with varying material strengths are head-on 
collisions with a relative velocity of 10 m/s. The basic relations 
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Fig. 8. Head-on collision with impact velocity 10 m/s. The ma- 
terial strength coefficients are So = 600 Pa and Tq - 600 Pa. The 
plot shows colour-scaled the density for each particle at simula- 
tion time f = 0.5 s after contact. 



equal sized ice grain agglomerates 




Fig. 10. Head-on collision with impact velocity 10 m/s. The ma- 
terial strength coefficients are So = 60 000 Pa and To - 6 000 Pa. 
The plot shows colour-scaled the density for each particle at sim- 
ulation time f = 0.5 s after contact. 



* 




Fig. 9. Head-on collision with impact velocity 10 m/s. The ma- 
terial strength coefficients are Eo = 6 000 Pa and To = 6 000 Pa. 
The plot shows colour-scaled the density for each particle at sim- 
ulation time t - 0.5 s after contact. 




Fig. 11. Head-on collision with collision velocity 10 m/s. The 
material strength coefficients are So - 60000 Pa and To = 
60 000 Pa. The plot shows a colour-scaled the density for each 
particle at simulation time t - 0.5 s after contact. 



for the compressive and tensile strengths and the bulk modulus 
remain unchanged, only the coefficients So and To are varied. 

The figureslTlfTOlshow colour-scaled density plots of the par- 
ticles at the end of each simulation. For the first of these simula- 
tions, we have used the same material properties as in sect. 15.11 
So = 600 Pa and Tq = 6 x 10^^ Pa (see fig. Q. Because the col- 
lision velocity is half the collision velocity from the simulation 
shown in fig. [1] the spheres do not experience the same dras- 
tic plastic deformation. In fact, they end up as two solid half 
spheres and no particles were pulled out of the spheres during 
the compression phase. Their relative velocity at the end are 
again vanishing. This picture changes drastically if the tensile 
strength coefficient To is reduced by one order of magnitude to 
600 Pa (see fig.[8j. Now the tensile forces are too weak to sustain 
the bodies during the contact phase. In the end, the two spheres 
stick together into a disc-like structure and many particles are 
ejected due to the fragmentation by plastic flow. If the compres- 
sive strength parameter So is increased to 6 x lO-' Pa = To instead 
(see fig.m, the material is strong enough to evade the fragmen- 
tation, only a few particles are thrown out, and the two spheres 
finally stick together, forming an elongated body. Interestingly, if 
the compressive strength is augmented even further to 6 x lO'* Pa 
(see fig. [Tol l, the bodies become too elastic and do not stick. 



The basic structure of the spheres is conserved, only in the im- 
pact contact region, the spheres are flattened. The objects rather 
bounce off' from each other, losing some particles. A nearly fully 
elastic rebound of the spheres happens for So = To = 6 x 10'* Pa 
(see fig.fTTTi. 

Again, the density of the spheres does not change signifi- 
cantly. Although their shapes change notably, they are not com- 
pacted effectively and their porosities stay mainly unchanged. 

6. Discussion 

T he simu l ation s presented in this work extend the investigation 
of lSironol (12004,) to collisions between equal-sized agglomerates. 
Sirono's main results regarding collisional growth and stickiness 
conditions can be summarised in the following way: Porous ag- 
glomerates can grow by collisions if the tensile strength is larger 
than the compressive strength, the shear strength is larger than 
the compressive strength and the Mach number is lower than 
0.04 for oblique impacts, which corresponds in our simulations 
to a maximum collision velocity of 3.1 m/s. Additionally, a dam- 
age restoration effect has to be included, otherwise the agglom- 
erates are totally fragmented. Since we do not use the damage 
model from Sirono, there is no need for a restoration effect. 
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Our new simulations with varying impact parameters indi- 
cate that collisions with b < 0.4 lead to drastic shape deforma- 
tions, the spheres finally form disc-like structures. Clearly, disc- 
like structures will withstand subsequent collisions even worse 
if the impact is perpendicular to the disc plane. On the other 
hand most (for head-on collisions nearly all) of the kinetic en- 
ergy is lost, which possibly leads to lower collision velocities for 
following encounters, thus favouring growth in subsequent colli- 
sions. The change in shape is lower for larger impact parameters 
b > 0.4, where less kinetic energy is lost. However, agglomer- 
ates emerging from oblique collisions tend to rotate and eventu- 
ally break if the tensile strength is not large enough. Then, plastic 
flow leads to fragmentation of agglomerates after the compres- 
sion phase. However, the tensile strength of a realistic porous 
protoplanetesimal might be higher in orders of magnitude than 
assumed here. 

One of the most striking simulation results is the miss- 
ing compaction of the agglomerates. Although they are rather 
porous objects (with 90 % porosity), their porosity does not 
change significantly during the impact. It is unclear, how the 
porosity in the agglomerates may be reduced and the bodies 
compacted. Probably, lower collision velocities and differing 
sizes of the colliding agglomerates lead to the incorporation of 
smaller bodies into large agglomerates, and a slight decrease of 
porosity. However, the velocity distribution of protoplanetesi- 
mals in the solar nebula at the Earth's orbital distance ranges up 
to several tens of m/s. One collision at such a high speed might 
destroy previous growth entirely. 

The varying material strength parameters which were inves- 
tigated in the second series of simulations show the strong influ- 
ence of the material properties on the collisional outcome. For 
the applied values we find the whole spectrum: sticking, no stick- 
ing but vanishing relative velocity, complete fragmentation, and 
elastic rebound. In contrast to Sirono's results, the only simula- 
tion that results in sticking uses Eo - To > Yq. Moreover, since 
the Mach number of our simulations is 0. 1 3 for the 1 m/s colli- 
sions and 0.26 for the 20 m/s collisions respectively, all simula- 
tions should lead to fragmentation according to Sirono's results. 

Altogether the results indicate that to prevent fragmenta- 
tion, a large compressive strength may compensate a low tensile 
strength (see fig.fTOli and vice versa (see fig.|2ll. 



7. Conclusion 

As long as we do not have more insight into realistic material 
properties of protoplanetesimals, it is cumbersome to give any 
accurate description about the formation of planetesimals by col- 
lisional growth. However, we have demonstrated that specific 
strength parameters can lead to growth and encourage the forma- 
tion of planetesimals even for collisions of equal sized objects, 
where previous investigations have found destruction. CuiTently, 
it is a major topic to establish realistic material parameters of 
protoplanete simals in laboratory experiments. Recent experi- 
mental data dBlum & Schrapleil2004t IWurm et al.ll2005allbl) pro- 
vide first insight in the properties of porous dust agglomerates 
which diff'er significantly from the properties applied in this 
study and may lead to further perception. We plan to recalibrate 
the existing SPH-model which was applied in this paper by per- 
forming comparison calculations of impact experiments. Thus, 
a realistic model for elastic and plastic behaviour of porous ag- 
glomerates will be developed. 

Additionally, our studies indicate that rotation is an impor- 
tant unexplored eff'ect which needs to be probed in more detail. A 



subsequent study will therefore focus on the collisions between 
rotating porous protoplanetesimals. 
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